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ABSTRACT 

Motivated by recent experimental results on charm physics we investigate implications of the 
updated constraints of new physics in rare charm meson decays. We first reconsider effects of the 
MSSM in c — > n7 constrained by the recent experimental evidence on Amo and find, that due 
to the dominance of long distance physics, D — > V-y decay rates cannot be modified by MSSM 
contributions. Then we consider effects of the extra heavy up vector-like quark models on the 
decay spectrum of — > 7r~^i~^i~ and — > K~^i~^i~ decays. We find a possibility for the tiny 
increase of the differential decay rate in the region of large dilepton mass. The i2-parity violating 
supersymmetric model can also modify short distance dynamics in c ^ u£~^£~ decays. We constrain 
relevant parameters using current upper bound on the — > TT~^i~^i~ decay rate and investigate 
impact of that constraint on the Df K'^i~^i~ differential decay dilepton distribution. Present 
bounds still allow small modification of the standard model differential decay rate distribution. 
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INTRODUCTION 



Recently the evidence for — oscillations has been reported by Belle and BaBar collab- 
orations [l[ 0]. Combining the measured quantities indicates nonzero AT as well as nonzero 
AniD 



x = Am^/rz) = 0.0087 ±0.003 . (1) 
These results immediately stimulated many studies (see eg. 

aaaaaa 

) . The obtained results 

for the relevant parameters describing — mixing are not in favor of new physics effects. 
However, they give additional constraints on physics beyond the standard model (SM) as already 
given in papers ^, ^. On the other hand, the study of rare D meson decays is not considered to be 
very informative in current searches of physics beyond the standard model IQ, U, 12, 13, 14, 1^, 1^, 
\y\ . as it is expected from "b" physics. Namely, most of the charm meson processes, where the flavor 
changing neutral currents (FCNC) effects might be present like c ^ u and cu cu transitions, are 



dominated by the standard model long-distance contributions 



Id, 



11 
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13 



14 



15|, 



13, 



Due to the GIM mechanism and smallness of the down-type quark masses the radiative c — > U7 
decay rate is strongly suppressed at the leading order in the SM [l^. The QCD effects enhance 



it up to the order of 10 ^ [19 
supersymmetric SM (MSSM) 



. New bounds on the mass insertion parameters within minimal 
1, [5] are derived using the — oscillations. We include into 



consideration the possible effect of MSSM with non-universal soft breaking terms on c ^ n7 along 



the lines of 



20 



211]. Although this approach leads to the enhancement of the SM value by a factor 



10, it is too small to give any observable effects in D ^ Vj decays {V is a light vector meson). 
The dominating long-distance (LD) contributions in the D — > V'j decays give the branching ratios 
of the order Br ~ 10~^ jiol. [l^. which makes the search for new physics effects impossible. 

Another possibility to search for the effects of new physics in the charm sector is offered in 
the studies of -D ^ Xi~^£~ decays which might be result of the c ui^£~ FCNC transition [3, 
Here X can be light vector meson V or pseudoscalar meson P. Within SM 
inclusion of renormalization group improved QCD corrections for the inclusive c — > ui^i~ gave 
an additional significant suppression leading to the rates T{c ue'^ e~) /T = 2.4 x 10"^*^ and 



r(c uh'^ii~)/Tdo = 0.5 X 10" 



-10 



22l | . These transitions are largely driven by a virtual photon at 



low dilepton mass mu = \J (p+ -|- while the total rate for D — > Xt^ I is dominated by the 
LD resonant contributions at dilepton masses m^^ = mp,muj,m^ 

New physics could possibly modify the dilepton mass distribution below p or distribution above 
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<j) resonance. In the case of -D — > Tri~^i~ there is a broad kinematical region of dilepton mass above 
(p resonance which presents an unique possibihty to study c — > ui'^i~ at high mu [l^ . 

The leading contribution to c ^ ui~^i~ in general MSSM with the conserved i?-parity comes 
from one-loop diagram with gluino and squarks in the loop 



HI, 



16 



23l ]. It proceeds via virtual pho- 



ton and enhances the c — > ui~^i spectrum at small mu. We find that bounds on the mass insertion 
parameters make the abovementioned enhancement in D rare decay 111, [3] to be negligible. 



Some models of new physics contain an extra up- 
FCNCs at tree level for the up-quark sector 



13 



25 



ike heavy quark singlet 



241 ] inducing the 



m, mm I, while the neutral current for the 



down-like quarks is the same as in the SM. The stringest bound on these models comes from the 
recent bound on Am in the — transition as given in 

HQ 

. In our calculation, we analyze 

how these bounds on the FCNC vertex cuZ affect the D — > Pi~^i~ decays. A particular version 
of the model with tree-level up-quark FCNC transitions is the Littlest Higgs model [29|. In this 



case the magnitude of the relevant c uZ coupling is even further constrained by the large scale 
f > 0{1 TeV) using the precision electroweak data. The smallness implies that the effect of this 
particular model on c — > ui~^£~ decay and relevant rare D decays is insignificant jisl . 

Among discussed models of new physics the super symmetric extension of the SM inclu ding the 
i?-parity violation is still not constrained as other new physics models. As noticed by [ll, 22 1 one 
can test some combinations of the i?-parity violating contributions in D"^ — > TT~^i^i~ decays. We 
place new constraints on the relevant parameters and search for the effects of new physics in the 
— > K^i~^i~ decays which might be interesting for the experimental studies. 



3d, 



3l|] and FERMILAB [32, \m collabora- 



There are intensive experimental efforts by CLEO 
tions to improve the upper limits on the rates for D — > Xi~^i~ decays. Two events in the channel 
^ .U*e- with m„ close .o have already been observed by CLEO fl. The other rare 
D meson decays are not so easily accessible by experimental searches, but with the plans to make 
more experimental studies in rare charm decays at CLEO-c, Tevatron and at charm physics sec- 
tions at present S-factories and in the future at LHC-b facilities makes the study of rare D decays 
more attractive. 

In order to compare effects of new physics and the standard model we have to determine size 
of the LD contributions. As in the case of — > Tr'^e~^e~ decay we use experimental data on the 
Dg nonleptonic decays accompanied by the vector meson dominance as we did in |l^ . 

The paper is organized as follows. In Section 2 we consider the impact of new charm mixing 
bounds on the c — > n7 decay. In Sec. 3 we study how new physics affects c — > ui~^i~ and 
D Pi'^i" decays. We present framework for calculating SD effects as well as the details of LD 
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calculations. In Sec. 4 we discuss our results and in Sec. 5 we make a short summary. 

II. c ^ M7 DECAY 

Given the recent observation of — mixing, we re-evaluate the possible effect of MSSM on 
c wy. Since the model with universal soft-breaking terms is known to have negligible effect {21]. 
we consider the model with non-universal soft breaking terms. We consider only the gluino exchange 
diagrams through {6^o)t,r.bt, mass insertions, since the remaining SUSY contributions can not have 



sizable effect 



2Q|, 



211 ]. The maximal value of {Si2)lr,rl insertion has been constrained by saturating 



X = Am£)/T = (4.8 it 2.8) x 10 with the gluino exhange in Q]. The results corresponding to the 
measured x = (8.7 ± 3) x 10~^ are shown in second column of Table HI Another constraint is 

obtained by requiring the minima of MSSM scalar potential do not break electric charge or color 
and that they are bounded from above {6i2)lr,b.l ^ V^mc/rriq Sj] , with values given in third 
column of Table [H The second constraint is obviously stronger for > 350 GeV, while Amu 
gives more stringent constraint for lighter squarks. Using {Si2)lr,rl < md 



GeV, m-c = 1.25 GeV and expressions from 



rriq, rriq = rUg 



2l[ | we get the upper bound 



r(c ^ uj)/rDO < 8 X 10^^ 



350 



(2) 



which is one order of magnitude larger than the standard model prediction 

r(c ^ uj)/Tuo = 2.5 X 10^8 m 



niq = rrig 


\{Si2)lr,b.l\ 


\{Si2)lr,rl\ 




from Amjj 


from stability bound 


350 GcV 


0.007 


0.006 


500 GcV 


0.01 


0.004 


1000 GeV 


0.02 


0.002 



Table 1: Upper bounds on mass insertions \{Si2)lr,rl\ from measured Am^) and stability bound 
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However, this possible SUSY enhancement by factor 10 would not affect the rate of the D V7 
decays, which are completely dominated by LD contributions with Br ~ 10~^ 



101, 



11 



The only window for probing the c — 
contributions are strongly suppressed 



7x7 enhancement remains the Be i?^7 decay, where LD 



35|. 
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III. 



d+l- AND D Pe+e- DECAYS 



SD Effects 



The c ui^l transition is driven by the low-energy effective Lagrangian 

i=7,9,10 

given in terms of four-quark operators 



(3) 



Qi = -^mcF^^ua^" {I + 75)c, 



<5io 



167r2 



167r2 



(4) 
(5) 
(6) 



and the corresponding Wilson coefficients Cy^g^io- F^^, is the electromagnetic field strength, while 
9l = |(1 1 are the left-handed quark fields. Wilson coefficients are taken at the scale /x = rric- 
Since we consider exclusive decay modes D — > Pi'^£~ we have to employ form factor description 
of the four-quark operators evaluated between two mesonic states. We use the standard parame- 
terization 



{P{k) I u^^il - 75)c I D{p)) = ip+ kYf+{q^) + {p- krUq^), (7) 
{P{k) I ua'^'^il ± 75)c I D{p)) = is(q^) [{p + kTq" - q^'ip + kf ± ie^''"^{p + kUfs] , (8) 

where P = vr"*" (K^) in the case of = (Df). The momentum transfer q = p — k is also the 
momentum of the lepton pair. For the /+ form factor we use the double pole parameterization of 
Ref. [SQ] 

/+(0) 



(1 — x){l — ax) 



(9) 



where x = q^/mj^t,, /+(0) = 0.617, and a = 0.579. We approximate s(g^) by f+{q'^)/niD, which 



is valid in the 
noticed in 38l . 



imit of heavy c-quark and zero recoil limit 



37l |. This relation can be modified as 



39[. However, in our case this modification cannot give significant effects. Finally, 



we arrive at the short distance amplitude for c ui~^£ decay 



A: 



SD 



-I —V^^V^b 



V2 



+ 



Cio 



u{p-)^v{p+) + --^u{p_)^-f^v{p. 



Uii). (10) 



2tt'^ m£) IGtt^ / IGvr^ 
p, p^ and p_ are the momenta of the initial D meson and the lepton pair in the final state, 
respectively. 
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Standard model 



The SM rate is dominated by the photon exchange, where c ^ n7 is a two loop diagram induced 

exchange Q, Q, Q, Q]. The effecti 



by effective weak vertex and a gkion 
is [l3] 



V:,VubCf = K:K.(0.007 + 0.0200(1 ± 0.2). 



tive Wilson coefficient 



(11) 



The remaining two Wilson coefficients are subdominant in the SM and we ignore them in further 
analysis. Cg is small due to the effects of the renormalization group, while the coefficient Cio is 
completely negligible in the SM [l^ . 



2. Models with extra heavy up vector-like quark singlet 



The class of models with an extra up-like quark singlet (EQS) naturally accommodate FCNCs 



at tree level 



13|, 



2M 



9 



cos 6w 



"^EM electromagnetic current, while the weak neutral current 



7/^ 



(12) 



(13) 



mixes up-type quarks 



29( 1 ■ where C/™" and D"^ are the quark mass eigenstates. The transition 
matrix to the mass eigenbasis for the up- type quarks is 4 x 4 unitary matrix T^, which causes 
tree- level FCNCs in the interaction term J/^gZ^ in the up sector. The mixing matrix contains only 



the elements of the last column of matrix 



n 



-@u&*T ^ 



(14) 



-QcOi i-ie^p -G^e^* - 
-etei -6*6* i-|etp -GtG^ 

\-6rG: -GtG* -GtGj* 1 - |QtIV 
The unitarity of the extended CKM matrix then implies that off-diagonal elements of 17 are non- 
zero, e.g. Que = -0«0c = ^udVcd + ^usVcs + ^ubVcb ^ 0. The low-energy effective description is 
encoded in Wilson coefficients Cg and Ciq. Relative to the negligible SM values, they are modified 
by the presence of an extra up-like quark: 



a 

Air . 
a 



0„c(4sin2 Ow-l) 



(15) 
(16) 
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The element Que of the up- type quark mixing matrix is constrained by the measurements of — 



mixing 



J] and using expression Am^ = 2 x 10 



n,.r < 2.8 X 10"^. 



(17) 



3. Minimal super symmetric SM 

The leading contribution to c — > ui~^i~ in general MSSM with conserved i?-parity comes from 
the gluino exchange diagram via virtual photon and significantly enhances c — > ui~^i~ at small 
m«. This MSSM enhancement can not be so drastic in hadronic decays, since gauge invariance 
imposes additional factor of m'j^ for D Pi'^i" decays, while D Vi~^i~ has large long distance 
contribution at small mu just like D — > V'j. 

In the MSSM with broken i?-parity (MSSM^), the c — > u£^i~ process is mediated by the 
tree- level exchange of down squarks [ll|. Integrating them out leads to the effective four-quark 
interaction 



eff 



E ^^(^L7'^Ci)(^"L7A). 



-3k 



are the CKM-rotated couplings between the L, Q and D supermultip^ 



ijk 

tial 



(18) 

ets in the superpoten- 



11[. In our notation the contribution to the Wilson coefficients is 

Vcb^ub^Cg = -v*f,VubSCiQ = — 5-!^ y ( — ^ I x'i2kKik: 




(19) 



where i = 1 (2) contributes to the e'^e (al!j^_) mode. The A'j^gfc ^^'^ -^ufe have been con 



strained from the charged current universality 
comes 



22 



4l|, while the strictest constraint on J2k '^22k^2ik 



22( 1 from the experimental limit Br(L'^ 7r'^/_f+/i ) = 8.8 x 10 



42]. We shall reanalyze 



the latter case, where LD physics generates the fair amount of experimental branching ratio. It 
is sensible to use an approach, where one takes into account the interference between LD and 
MSSM^ part of the amplitude to constrain the couplings of the MSSM^. 

B. Long distance contributions in — > P£~^£~ 



Knowledge of the LD contributions is crucial, if we want to isolate short distance physics in 
the decays of type D — > P£~^£~ . Following procedure described in 1^ we consider long distance 
contributions by employing the resonant decay modes, in which D first decays to P and a virtual 
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neutral vector meson Vq, followed by decay of Vq — > 7 — > t^^ ■ First stage of the decay is controlled 
by effective weak non-leptonic Lagrangian 

G p 

q=d,s 

The effective Wilson coefficients on the scale are 



4° = y^^y*^ [ain7'^(l - ^^)qq^^{l - 75)0 + a2«7''(l - l^)cq^^{l - 75)^7] (20) 



22] 



ai = 1.26, 02 = -0.49. 



(21) 



The flavour structure of (j20p allows Vq to be either p, uj or (j). Since branching ratios of separate 
stages in the cascade are well measured, we will not work in a particular theoretical model, but 
will instead try to make the best use of experimental data currently available. Here we follow the 



lines of Ref. 

dr 



13t |. For a cascade, we write [43 [ 



1. 



Here T rj^^KVoiff) and Ty^-^^i^+n- denote decay rates if Vq had a mass and these rates are 



,2 r2 ^ Vo^l+t 



-iq'). (22) 



known experimentally only at = mvo- Since the resonances Vq = p,uj,(j) are relatively narrow 
(Fyg <C rny^), the following relation approximately holds 



Br [D PVo Prt] = Br [D PVo] x Br [Vq . 

The phenomenological amplitude ansatz that reproduces the above behaviour is then [1 

ava 



(23) 



-(p„)^+(p+)] =e^^-o 



(24) 

The only assumption we made here is that coefficient avo is independent of q^. We included the 
phase (pVo explicitly, so that ayg is real and positive number. 

1. D+ TT+£+i- 



For the right side of Eq. (j23p we use experimental data (Table ln|). and use ansatz (|24p to extract 
unknown parameters avo- o-p = 2.94 x 10~^, = 4.31 x 10~^. Decay mode — > vr+u; has not 



mode 


D+ -f 71+ p 


D+ TT+UJ 


D+ n+cj) 


BrxlO^ 


1.07±0.11 


< 0.34 


6.50 ±0.70 



Tabic II: Branching ratios of decays of D+ meson to the intermediate resonant states 



44| 
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been measured yet, but we can relate a^o and its phase to the well-measured contribution of the p 
resonance assuming vector meson dominance as in jisl . Relative phases and magnitudes of the reso- 
nances are extracted by considering the decay mechanism, controlled by the weak Lagrangian (j20p 
and electromagnetic coupling of Vq to photon. Then flavour structure of the resonances deter- 
mines relative sizes and phases of resonant amplitudes. Detailed analysis has already been done 
in Ref. The relative phases of p and lo contributions are found to be opposite in sign, while 
for the ratio of they magnitudes it was found that a^/ap = 1/3. Also the phases of p and (j) are 
opposite. Thus the final LD amplitude (up to the phase) becomes 



iM 



LD 



3 g2 _ 7^^2 _|_ im^Y^ J -m\ + im^T^ 



(25) 



2. D+ K+i-£+ 



Experimental data is not as rich as in the case of non-strange charmed meson decays (Table 
Contributions of p and uj are related like in the case of meson, namely a^jap = 1/3 with 



mode 


D+ -> K+p 




K+uj 


D+ K+(j) 


BrxlO^ 


2.60 ± .70 




< 0.50 



Table III: Branching ratios of Df meson to the intermediate resonant state 



44|. 



opposite relative phase between them. In the same way as for the non-strange decays, we determine 
Op = 6.97 X 10~^. However, the contribution of the (j) resonance is only limited from above 
by experimental data and we have to rely on a theoretical model. Consequently, the total LD 
amplitude is a sum of two terms: 



A 



LD 



1 



1 



g2 — + iTTT-oT^ 



(26) 



We calculate the <j) part of (I26p using the vector meson dominance (VMD) assumption, where the 
intermediate (p contributes by decaying into a virtual photon, which further decays to the lepton 
pair. Both oi and a2 parts of the non-leptonic Lagrangian (j20|) can generate the flavour quantum 
numbers of (j) and K~^. The oi part connects initial state to (j) through a charged current 
{sc)v-A, while the {us)v-a creates the K'^ out of vacuum. Neutral currents (the 02 part) do the 
opposite: Df and — > (/>. Utilizing the Feynman rules of (j20p . and the VMD hypothesis 
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we arrive at the (f) contribution to the LD amphtude 

Q [Q ~ i^ff, ~r iiT^4>^ 4>) 

(27) 

However, the phase of relative to the rest of the amphtude (j26p remains to be free. The 



P — > y transition ^ is described by the form factor Aq, which we take from [45|. In our 
calculations we consider gauge invariant amplitude for Df K~^£~^£~ in which dependence 
is cancelled. 



IV. RESULTS 



Using the approach, described in Sec. 3, we analyze impact of short distance physics on long- 
distance resonant background. Since the SD contribution of SM is completely overshadowed by 
LD, we will only consider the EQS and MSSM^ models of new physics, to see if the experimental 
searches for them are still viable in — > Xi~^i~ decays. Current constraints on EQS model coming 
from the — mixing already indicate the dominance of LD contributions in the total decay 
rate. On the other hand, the contribution of MSSM^ is not as constrained and one should still 
see the deviations from the LD contribution away from the resonant region of the phase space. 

We shall analyze the dilepton squared mass {q'^ = rn'j^) distribution of the branching ratios. We 
fix free phases in the amplitude in a way which maximizes branching ratio for the considered decay 
mode. 



A. D+ n+e+i- 



mode 


LD 


extra heavy quark 


LD+extra heavy quark 


MSSM^ 


LD+MSSM^ 


D+ TT+e+e.- 


2.0 X 10"*^ 


1.3 X 10-^ 


2.0 X 10"^ 


2.1 X 10-^ 


2.3 X 10-*^ 


D+ TT+fl+fJ.- 


2.0 X 10"^ 


1.6 X 10-9 


2.0 X 10"^ 


6.5 X 10"^ 


8.8 X 10"^ 



Table IV: Total branching fractions of the L)+ 7r+£+£^ modes. In the first column (LD) are only long- 
distance BRs. The remaining four columns give maximal contributions of the SD physics models alone and 
also combined contributions of the SD and LD physics. 



Branching fractions are listed in Table IIVI Clearly, the EQS model contribution is too small to 
be observed (Fig. [1]). 
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LD+extra quark sing. 

LD 

extra quark sing. 




q^ (GeV^) 



LD+extra quark sing. 

LD 

extra quark sing. 




q^ (GeV^) 



Figure 1 : Distributions of the maximal branching ratios in the model with extra quark singlet for the decay 
modes — + 7r+e+e^ (left) and 7r+/i+/i^ (right). Full line represents the combined LD and SD 

contributions. 



On the other hand, the MSSM.^ gives a shght increase to the mode with electrons. Deviation 
from the LD amphtude is pronounced in the region without resonances, where mu < rup or 
mu > fn^f) (Fig. [21 left). However, the most promising mode is the channel with muons. The long- 



5 




q^ (GeV^) 



0.001 
1e-04 
1e-05 
1e-06 
1e-07 
1e-08 
1e-09 
1e-10 



0.5 





LD+MSSM-R 




LD 




MSSIVI-R 




ti^i— ^^^^ 



1.5 

q^ (GeV^) 



2.5 



Figure 2: Distributions of the maximal branching ratios in the MSSM_^ model for the decay modes D'^ 
Tr~^e'^e~ (left) and 7r+/i+/i^ (right). Full line represents the combined LD and SD contributions, and 

it corresponds to the experimental upper bound BR(D^ fi'^ ii' ") = 8.8 X 10"^ on the right plot. 



i3 (8.8 X 10-6) 



distance contribution (2.2 x 10" ) is a fair share of the experimental upper bound 
and should be taken into account together with the short distance part, when one is constraining 
the Wilson coefficients. The difference is not big, i.e. when we drop the LD part we get for the 
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bound I V^.V^feC'g j^qI < 27, while the analysis with LD part included gives 

iKlKfeCg^.iol < 23. (28) 

The latter bound is a maximum with respect to the free relative phase between the LD and 
MSSM^ parts of the amplitude. Although the inclusion of the LD term does not make substantial 
difference, it will grow rapidly as experimental bound is approaching 2.2 x 10~^. All the branch- 
ing ratios concerning MSSM^ and muons in the final state (Table IIV|VP and their kinematical 
distributions (Fig. EH]) use the bound ([28]) . 

B. D+ ^ K+i+e- 

The branching ratios contributions are summarized in Table jVj Again, the EQS model has 
negligible effect (Fig. [3]). MSSM^ has a notable effect, especially in the /^'^/U" mode, where it 
increases branching ratio by an order of magnitude (Fig. H]). In this case, the MSSM^ overshadows 
the LD contribution throughout the phase space, except in the close vicinity of the LD resonant 
peaks. 



mode 


LD 


EQS 


LD+EQS 


MSSM/? 


LD+MSSM^ 


D+ K+e+er 


6.0 X 10-^ 


5.4 X 10-1° 


6.0 X lO-'^ 


9 X 10-** 


7.6 X lO"'^ 


D+ ^ K+^,+^l- 


6.0 X lO"'^ 


6.2 X 10-1° 


6.0 X lO-'^ 


2.6 X 10-s 


3.6 X 10-6 



Table V: Total branching fractions of the Z?+ K^£^£^ modes. In the first colmim (LD) are only long- 
distance BRs. The remaining four columns give maximal contributions of the SD physics models alone and 
also combined contributions of the SD and LD physics. 

V. SUMMARY 

Recently observed — mass difference constrains the value of tree-level flavor changing 
neutral coupling c — > uZ, which is present in the models with an additional singlet up-like quark. 
We have studied the impact of this coupling on rare 'k'^1~^£~ and Df K~^£~^l~ decays, 

where its effects are accompanied by the long distance contributions. We have determined long- 
distance contributions in Df — > K'^l~^i~ following the same phenomenologically inspired model 
as it has been done previously in the case of 7r~^i~^£~ . We find that the effect of new extra 

singlet up-like quark is too small to be seen in dilepton mass distributions for both decay modes. 
In our previous study we have considered forward-backward asymmetry in the — > p^i~^i~ and 
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LD+extra quark sing. 

LD 

extra quark sing. 




q^ (GeV^) 



2.5 



LD+extra quark sing. 

LD 

extra quark sing. 




q^ (GeV^) 



2.5 



Figure 3: Distributions of the maximal branching ratios in the model with extra quark singlet for the decay 
modes Df /\+e+e^ (left) and Df K^iJ^n^ (right). Full line represents the combined LD and SD 
contributions. 




2.5 



0.001 



1e-04 



1e-05 



1e-06 



1e-07 



1e-08 





LD+MSSM-R 

LD 

MSSM-R 







0.5 



1 1.5 

q^ (GeV^) 



2.5 



Figure 4: Distributions of the maximal branching ratios in the MSSM^ model for the decay modes Df 
K^e^e~ (left) and Z?+ ji^ ji'^ (right). Full line represents combined LD and SD contributions. 



found very small effect. New constraint reduces that asymmetry even more, making it insignificant 
for the experimental searches. 

Present constraints on mass insertions in MSSM with conserved i?-parity still allow for increase 
of c — > ii7 rate by one order of magnitude. For the same reason MSSM could significantly increase 
c — > ui~^i~ rate at small mu. However, this MSSM enhancement is not drastic in D decays, since 
D and D Vi~^i~ have large long distance contributions for small rriu, while D Pi'^£~ 

rate is multiplied by factor of mj^ due to gauge invariance. 

The remaining possibility to search for new physics in rare D decays is offered by the MSSM 
models which contain R-parity violating terms. We reinvestigate bounds on the combinations of 
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these parameters in Tr~^£^£ by including the long-distance effects. Using current upper 

bound on the rate for — > -rr^l'^i"' we derive new bound: 

E(^^^)'a'2U.A^2.< 0.0041 (29) 

fc = l V '^R / 

Since at Tevatron there are plans to investigate — > K^i~^i~ decay we use upper bound (j29p 
and calculate dilepton invariant mass distribution. This bound still gives small increase of the 
dilepton invariant mass distribution for the larger invariant dilepton mass, making it attractive for 
the planned experimental studies. 
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